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摘 要 :为 了 预测 三 维 编织 C/C 复合 材料 的 谊 曲 失效 行为 ,基于 多 尺度 渐进 展开 理论 ,结合 细 观 渐 
进 损伤 模型 ,建立 了 三 维 编织 C/C 复合 材料 宏 细 观 多 尺度 分 析 模 型 。 通 过 商业 有 限 元 软件 
ABAQUS 用 户 子 程序 UMAT 的 二 次 开发 ,在 宏观 结构 有 限 元 分 析 中 实时 调用 细 观 单 胞 模型 进行 细 
观 渐进 损伤 分 析 , 实现 了 宾 细 观 尺 度 之 间 交 互 式 信息 传递 和 多 尺度 损伤 模拟 。 利 用 上 述 模型 对 三 
点 弯曲 载荷 下 三 维 编织 C/C. 复合 材料 梁 的 渐进 损伤 和 失效 过 程 进行 了 模拟 ,预测 了 梁 的 载荷 - 挠 度 
) 曲 线 和 弯曲 强度 ,并 与 实验 结果 进行 了 对 比分 析 , 验 证 了 基于 多 尺度 方法 的 三 维 编织 C/C 复合 材料 
弯曲 强度 预测 模型 的 有 效 性 ,为 此 类 材料 及 结构 失效 分 析 提 供 了 一 种 手段 。 
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Bending failure characteristics of 3D braided C/C composites 
based on multi-scale asymptotic method 
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Abstract :In order to predict the bending failure behavior of the three-dimensional (3D) braided C/C 
composites,a multi-scale model was established based on the multi-scale asymptotic method. The repre- 
sentative volume element model for the material was called at each material integration point of the macro 
structure , the meso-stress fields for fiber yarns and matrix were computed , and the progressive damage a- 
nalysis was conducted at meso level. The information exchanges between different scales and progressive 
damage were accomplished through a user developed subroutine UMAT in ABAQUS. The load-deflection 
curve , damage propagation , and bending strength of the 3D braided C/C composite beam under three-point 
bending were predicted by the developed model. The predicted results were in good agreement with the ex- 
periment results , demonstrating the validity of the developed model, which provides an effective tool for 
failure analysis of such kind of composite structure. 
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三 维 编织 C/C 复合 材料 除了 具有 传统 C/C 复 
合 材料 高 比 强 .高 比 模 以 及 耐 高 温 等 优点 外 ,还 具有 
抗 分 层 . 抗 热 震 以 及 性 能 可 设计 等 优点 ,在 航空 航天 
等 高 科技 领域 得 到 广泛 应 用 ,常常 作为 大 型 固体 发 
动机 喷 管 喉 部 ,高速 飞行 器 到 前 缘 以 及 导弹 端 头 帆 
等 上 0 。 作 为 新 型 高 温 热 结构 材料 ,三维 编 织 C/C 复 
合 材料 制备 工艺 复杂 、 周 期 长 ,开展 性 能 预报 研究 对 
于 三 维 编织 C/C 复合 材料 结构 设计 和 材料 开发 具 
有 重要 意义 。 

目前 ,已 经 有 很 多 学 者 对 三 维 编织 复合 材料 的 
拉 伸 . 剪 切 以 及 黏 弹性 性 能 预报 进行 了 研究 5 。 在 
强度 预报 方面 ,Aubard 等 "建立 了 宏观 弹 塑 性 损 
伤 模型 预报 了 三 维 四 向 编织 C/C. 复合 材料 的 单 轴 
拉 伸 非 线性 力学 行为 。 文 献 [11-12] 利 用 连续 介质 
损伤 理论 研究 了 三 维 细 编 穿刺 C/C 复合 材料 在 复 
大 强 荷 下 的 非 线性 应 力 应 变 关系 。 文 献 [ 13-15 ] 通 
二 对 纤维 束 编织 结构 进行 适当 简化 ,采用 统计 理论 
VER T HERR C/C 复合 材料 的 拉 伸 强度 。 三 维 
编织 C/C 复合 材料 既是 一 种 材料 ,也 是 一 种 结构 ， 
传 迪 的 宏观 强度 理论 将 三 维 编织 复合 材料 等 效 为 均 
命 他 材料 ,难以 反映 材料 细 观 编织 结构 。 近 年 来 ,一 
ney 
梅 现 细 观 模型 ,采用 有 限 元 方法 研究 三 维 编织 C/C 
复 容 材料 力学 性 能 。 文 献 [16-17 ] 利 用 周期 性 单 胞 
模型 ,对 纤维 束 和 基体 分 别 采 用 细 观 强度 准则 ,分 别 
预报 了 细 编 穿刺 C/AC 复合 材料 和 三 维 针 刺 CIC 复 
合 靶 料 的 宏观 拉 伸 强度 ,分 析 了 纤维 束 /基体 界面 的 
影响 。 文 献 [18-19] 采 用 uCT 三 维 重 构 技术 建立 了 
REZAR C/C 复合 材料 单 胞 模型 ,利用 渐进 损 
伤 方法 研究 了 材料 制备 孔 阶 缺 陷 和 含量 对 三 维 正 交 
编织 C/C 复合 材料 单 轴 拉 伸 失 效 行为 的 影响 。 陈 
波 等 2 考虑 了 纤维 束 在 高 温 氧化 环境 下 的 性 能 退 
化 ,利用 渐进 损伤 方法 建立 了 高 温 氧 化 环境 下 三 维 
四 向 编织 C/C 复合 材料 拉 伸 强度 预报 模型 ,研究 了 
氧化 环境 对 材料 拉 伸 强度 和 损伤 扩展 的 影响 。 上 述 
研究 利用 单 胞 模型 从 细 观 角度 研究 了 三 维 编织 C/C 
复合 材料 单 轴 拉 伸 或 压缩 失效 行为 ,对 于 在 弯曲 载 
荷 下 三 维 编织 C/C. 复合 材料 的 失效 行为 的 预报 研 
究 较 少 。 在 弯曲 载荷 作用 下 ,三 维 编织 C/C 复合 材 
料 结构 的 宏观 应 力 应 变 旦 非 均匀 分 布 ,利用 单 胞 模 
型 难以 施加 载荷 和 边界 条 件 。 如 果 对 宏观 结构 全 域 
进行 细 观 结构 建 模 和 分 析 宏 观 与 细 观 之 间 的 尺度 差 
异 ,又 会 带 来 数值 模型 庞大 .甚至 难以 计算 等 问题 。 
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针对 上 述 问题 ,本 研究 基于 渐进 展开 多 尺度 理 
论 ,建立 了 三 维 编织 C/C 复合 材料 结构 失效 的 宏 细 
观 多 尺度 分 析 方 法 。 对 于 宏观 结构 采用 均 质 化 材料 
模型 ,在 宏观 积分 点 处 采用 能 够 反映 细 观 编织 结构 
的 单 胞 模型 ; 对 于 纤维 束 和 基体 采用 细 观 渐进 损伤 
模型 ,利用 商业 有 限 元 软件 ABAQUS 材料 子 程序 
UMAT 二 次 开发 ,实现 了 宏 细 观 之 间 交 互 式 信息 传 
递 和 多 尺度 损伤 模拟 ;预测 了 三 点 弯曲 载荷 下 三 维 
编织 C/C 复合 材料 梁 的 渐进 损伤 和 失效 行为 ,并 与 
实验 结果 进行 了 对 比 验 证 。 


1 方法 和 模型 


1.1 多 尺度 渐进 展开 方法 


如 图 1 所 示 , 对 于 具有 周期 性 微 结构 的 三 维 编 
织 C/C 复合 材料 ,定义 材料 宏观 域 为 2, 宏观 尺度 坐 
标 系 为 x; 定 义 其 细 观 周期 性 单 胞 域 为 Y, 细 观 尺度 
坐标 为 y, 宏观 和 细 观 太 度 之 比 为 = =y/x。 图 中 宏 
观 结构 坐标 系 * 方 向 对 应 细 观 单 胞 模型 的 六 方向 ， 
宏观 ,方向 对 应 细 观 单 胞 模型 的 y 方 向 。 
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图 1 三 维 编织 C/C 复合 材料 宏 细 观 结构 及 其 
周期 性 单 胞 模型 


Fig.1 Macro-micro structure and periodic unit cell 
model of three-dimensional braided C/C composite 
将 结构 的 宏观 位 移 wu 展开 成 关于 2 种 尺度 之 比 
e 的 渐 近 展开 式 为 
u*(x) zu(x,y) =u (x,y) * eu? (x,y) « & 
(1) 
同 理 , 可 将 应 变 展开 成 关于 e 的 渐进 展开 式 为 


1 (- 
ej(u’) = D (xy) +e (x,y) teet) (x,y) te 


(2) 
根据 几何 方程 .本 构 方 程 以 及 弹性 问题 平衡 方 
程 , 令 关于 e 老 的 各 阶 系数 为 0, 得 到 一 系列 摄 动 方 


程 , 即 
dc^? 
eg, Ma (3a) 
9y; 
E 0m; 0n; E 
ME. 0x; ðY; a jam) 
oc e „20 7 
0 
€ : m * 35 +f: = 0 (3c) 
同时 ,得 到 关于 宏观 位 移 (uW ) 和 一 阶 细 观 位 移 
(un) 的 控制 方程 为 
(0) a 
a(i m ) =0 (4) 
ay; "V Oy, ðN 


式 (4) 联 系 着 宏观 位 移 和 一 阶 细 观 位 移 ,并 且 可 以 
证 明 一 阶 细 观 位 移 具 有 了 周期 性 , 且 其 解 可 以 表 
示 为 
NR E " oui 

U; 二 e (5) 
ip XL 为 宏观 位 移 与 细 观 位 移 之 间 的 关联 矩阵 ， 
s 作 单 胞 特征 位 移 场 。 

习 根 据 式 (4) 、 式 (5 ) 及 弹性 本 构 关 系 ,可 以 得 到 
aei uote FEME Dx: 和 单 胞 一 阶 细 观 应 力 oy 
的 控制 方程 , 即 


3 kl (0) 
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< 可 以 看 出 ,结构 内 部 某 一 点 的 细 观 应 力 可 以 分 
AO WAY :一 部 分 是 由 该 点 处 组 分 材料 的 自身 刚度 
引起 ; 另 一 部 是 由 于 材料 细 观 结构 引起 的 附加 应 力 。 
HARG) MRG) 导 到 求解 均匀 化 宏观 位 移 场 
u O 的 控制 方程 为 

o) 
EHE] m 
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从 均匀 化 系数 的 定义 可 以 看 出 ,均匀 化 系数 即 
宏观 结构 的 等 效 刚度 ,在 得 到 宏观 结构 的 等 效 刚度 
后 ,就 可 以 求解 宏观 结构 的 位 移 场 。 可 以 看 出 ,求解 
9 匀 化 等 效 刚度 系数 最 关键 的 是 求解 单 胞 特征 位 移 
Mi。 这 可 以 采用 有 限 元 法 ”进行 数值 求解 。 


1.2 细 观 损伤 模型 
通过 单 胞 细 观 控制 方程 式 (7) 可 以 获得 宏观 结 


Xi Jar (9) 
ðY, 


构 所 对 应 的 细 观 应 力 ,结合 一 
可 以 建立 细 观 损伤 模型 。 

三 维 编织 C/C 复合 材 的 细 观 组 分 材料 主要 是 
纤维 束 和 基体 。 假 设 纤 维 东 为 正 交 各 向 异性 单 向 复 
合 材料 ,在 弯曲 载荷 作用 下 ,纤维 束 的 主要 破坏 模式 
有 :纤维 束 纵向 拉 伸 破坏 、 纤 维 束 纵向 压缩 破坏 、 纤 
维 束 横向 拉 伸 破坏 和 纤维 束 横向 压缩 破坏 。 为 简单 
起 见 ,采用 最 大 应 力 准 则 作为 上 述 4 种 破坏 模式 的 
损伤 准则 ,同时 取 拉 伸 和 压缩 破坏 准则 中 的 较 大 值 
作为 不 同方 向 的 损伤 判 据 。 


定 的 细 观 强度 准则 就 
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式 中 :7 TZ 分 别 表示 直角 坐标 系 下 纤维 束 的 纵向 
CL) 和 横向 (了 ,2Z) o (1I S LT, Z) 4E IEEE 2E R3 
AFERE; E Fi S L,T,Z) TIAA ER 3 
个 方向 的 拉 伸 和 压缩 强度 。 

对 于 基体 材料 , 假设 为 各 向 同性 材料 ,采用 式 
(11) 所 示 Von-Mises 等 效应 力 强度 准则 作为 初始 损 
伤 判 据 。 


alo -0,) +(o -os) +(01-0;) ] 


=1 (11) 
式 中 :4 为 基体 初始 损伤 方程 值 ;o, o. o. 为 基体 
材料 的 3 PEJ; E, 为 基体 材料 的 拉 伸 强度 。 

纤维 束 和 基体 组 分 材料 中 应 力 一 旦 满足 上 述 失 
效 准 则 后 , 则 意味 着 组 分 材料 发 生 损 伤 。 在 有 限 元 
模拟 中 ,通常 采用 折 减 其 刚度 的 方法 模拟 损伤 。 根 
据 刚 度 折 减 因 子 的 不 同 ,主要 有 固定 的 折 减 因子 和 
化 的 折 减 因子 两 种 折 减 方案 。 本 研究 采用 与 

组 分 材料 断裂 能 .单元 特征 尺度 和 局 部 应 力 相 关 的 
伤 演化 方程 ,引入 3 个 损伤 变量 d, dy 和 
dz ,描述 纤维 束 的 渐进 损伤 过 程 。 


dl -pop -Cue A - DE] an 


1 
d -z| E T0 og] (13) 


- Gs ( hs 1t] (14) 
式 中 :dj,dr 和 dz 分 别 为 纤维 束 纵 向 (ZL) 和 2 个 横 


1 
d, =1 -一 ex 
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HCT Z) 损伤 变量 ;由 ,rr 和 由 z 分 别 为 纤维 束 3 个 
方向 初始 损伤 方程 代数 值 ;Ci Co IE Cs 分别 为 纤维 
RI 个 方向 的 正 轴 刚 度 ;eny、ewy 和 ss, 分 别 为 纤维 
3&3 个 方向 的 破坏 应 变 ;LL 为 单元 特征 尺度 , 取 单元 
体积 的 3 次 方 根 ;Cr 和 Gn 分 别 为 纤维 束 和 基体 断 


eu 
裂 能 。 


对 于 基体 损伤 ,损伤 变量 d,, 通过 以 下 方程 获得 
1 L, 
d,-1- 5| -Cuens( Pn- 1 | (15) 


WP id, 为 基体 损伤 变量 ;es /为 基体 破坏 应 变 。 
假设 损伤 主 向 和 材料 主轴 重合 ,将 损伤 变量 d, 
(=Z,7,Z 7m) 引 入 组 分 材料 刚度 矩阵 C ,使 材料 刚 
度 随 着 损伤 的 发 展 而 逐渐 衰减 ,对 于 纤维 束 而 言 ,这 
个 含 损伤 变量 的 刚度 矩阵 C (d) 可 以 表示 为 
T Phi. 25C bbCs 0 0 0 
bCa dus. 0 0 0 
bCa 0 0 0 
0 
0 


sym brCa 0 


bzCss 


(16) 
xe, -1 -di;0 =1 -dr;bz =1 -dz;sym 表示 对 
FORE 


L3. 多 尺度 有 限 元 模型 


P EE 上述 突现 和 细 观 多 尺度 分 析 方法 ， 采用 
ABAQUS 商业 有 限 元 软件 建立 如 图 2 所 示 的 三 维 纺 
AX C/C 复合 材料 宏观 和 细 观 多 尺度 有 限 元 模型 ,分 
析 名 型 包括 细 观 有 限 元 模型 和 宏观 有 限 元 模型 。 宏 
观 模 型 为 长 x 宽 x 高 =60 mm x 14 mm x6. 8 mm 的 
三 维 编 织 C/C 复合 材料 梁 , 跨 距 为 50 mm; 采 用 8 35 
点 六 面体 单元 C3D8 进行 有 限 元 网 格 划 分 ,在 梁 中 
部 支撑 点 等 高 应 力 区 域 进行 局 部 网 格 细 化 ;边界 条 
件 为 深 的 一 端 固 文 ,一 端 简 支 ,通过 面 接触 将 载荷 施 
加 到 梁 上 表面 中 间 部 位 的 网 格 节点 上 。 

根据 三 维 编织 C/C. 复合 材料 的 周期 性 细 观 结 
构 , 在 三 点 弯曲 梁 的 宏观 有 限 元 模型 的 每 个 单元 积 
分 点 处 ,提取 能 够 反映 三 维 编织 C/C. 复合 材料 细 观 
编织 结构 的 最 小 重复 性 单 胞 。 单 胞 模型 的 长 x S x 
高 =5.4 mm x3.2 mm x3.3 mm。 单 胞 模型 包括 纤 
维 束 增强 相 和 基体 相 。 采 用 8 节点 六 面体 C3D8 单 
元 进行 网 格 划分 , 共 划 分 72 420 个 单元 ;纤维 束 单元 
和 基体 单元 为 共 节 点 连接 。 关 于 单 胞 模型 的 具体 建 
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图 2 三 维 编织 C/C 复合 材料 梁 宏 细 观 多 尺度 模型 
Fig.2  Macro-meso scale models of three-dimensional 


braided C/C composite beam 
1.4 周期 性 边界 条 件 


在 利用 细 观 问题 控制 方程 求解 单 胞 特征 函数 x 
时 , 单 胞 需要 满足 式 (17) 所 示 周 期 性 边界 条 件 。 假 
设 单 胞 的 形状 为 六 面体 , 式 (17a) 表示 单 胞 每 两 个 对 
面 中 位 置 相同 点 处 的 特征 位 移 相 等 ; 式 (17b) 表示 
单 胞 8 AA ABCDEF G 和 互 处 的 特征 位 移 
为 0。Xia 等 ”给 出 这 种 特殊 形式 的 边界 条 件 的 实 
现 方法 。 本 研究 利用 ABAQUS 商业 有 限 元 软件 中 
多 点 耦合 位 移 约 束 方 程 ( Multi-point coupling) 实现 
上 述 周期 性 边界 条 件 。 


xi (y+Y) on ðY (17a) 

xi (y) =0 on gy" (17b) 
RP ac O) 为 单 胞 位 移 ;y 为 细 观 坐标 ;了 为 单 胞 
域 ;97 为 单 胞 边界 ;97™ 为 单 胞 顶点 。 


1.5 组 分 材料 性 能 


三 维 编织 C/C 复合 材料 采用 T300 »x T 2H RIS 
青 碳 基体 制备 而 成 ,纤维 丝 和 碳 基 体 的 性 能 参数 见 
表 1。 相 关 性 能 参数 取 自 文献 [24-27 ] ,纤维 束 纵 癌 
压缩 强度 取 拉 伸 强 度 的 60%'”。 需 要 指出 的 是 ,经 
过 2 500 % 高 温 热处理 后 ,碳纤维 模 量 有 很 大 提高 ， 
强度 则 有 一 定 的 下 降 。 纤 维 束 的 性 能 参数 可 以 根据 
纤维 丝 .基体 性 能 及 纤维 体积 分 数 利 用 Chamis 理论 
公式 ”| 复合 得 到 。 三 维 编织 C/C 复合 材料 轴 向 纤 
维 束 体积 分 数 为 80% , 径 向 纤维 束 体 积分 数 为 
60% ,纤维 束 的 性 能 参数 如 表 2 所 示 。 
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RI 纤维 丝 和 基体 性 能 
Tab.1 Fiber yarn and matrix properties 

项 目 En / GPa Ep / GPa x Gy / GPa Es /GPa F,/MPa F,/MPa 

碳纤维 丝 350 40 0.26 24 14.3 1 600 960 

碳 基体 12 = 0.33 - = 35 60 

表 2 纤维 束 性 能 参数 
Tab.2 Fiber bundle performance parameters 
AR E,/GPa ~ E,/GPa ni Go/GPa ^ Gy4/GPa F! /MPa F/MPa Fi/MPa ~ Fs/MPa 

轴 向 纤维 束 282.4 32.09 0.274 16.48 11.45 1 280 768 32.68 56.03 
径 向 纤维 束 214.8 26.21 0.288 12.16 9.60 960 576 30.72 52.67 


2 结果 与 讨论 


| 


宏 细 观 应 力 分 析 


Cl 3 给 出 了 三 维 编织 C/C 复合 材料 三 点 弯曲 
S33 应 力 分 布 情况 (加 载 位 移 0. 14 mm ,载荷 为 
6052N) 。 $33 应 力 是 指 6 个 应 力 分 量 中 沿 着 图 中 从 
EA Z 方向 的 应 力 。 由 图 3 可 见 ; 在 弯曲 载荷 作用 
edet rds ton) Jets Ed 
ZOI S33 应 力 分 量 为 负 值 ,下 表面 S33. 应力 分 量 为 
, 梁 中 性 面 S33 应 力 近似 为 0; 宏观 应 力 分 布 答 

点 弯曲 梁 的 一 般 变 形 规律 。 
Ea 给 出 了 三 维 编织 C/C 复合 材料 三 点 弯曲 
奖 外 下 表面 应 力 最 大 点 处 对 应 单 史 的 833 细 观 应 力 


S.S33 

(Avg: 7596) 
-1.878e+07 
-3.662e+07 


(a) 上 表面 


2.2 宏观 载荷 - 挠 度 曲线 

图 5 为 预测 的 三 维 编织 C/C 复合 材料 三 点 弯 
曲 梁 的 载荷 - 搁 度 曲线 和 实验 测试 曲线 。 由 图 5 可 
见 :基于 多 尺度 模型 的 预测 曲线 和 实验 曲线 总 体 趋 


8,833 
(Avg: 7596) 


场 。 由 图 4 可 见 , 三 点 弯曲 载荷 下 梁 上 表面 对 应 单 
胞 的 S33 细 观 应 力 值 为 负 , 下 表面 对 应 单 胞 的 S33 
细 观 应 力 值 为 正 。 这 是 由 于 弯曲 载荷 下 梁 的 上 表面 
受 压 ,下 表面 受 拉 引起 的 ;同时 , 轴 向 纤维 束 中 的 应 
力 水 平 远 高 于 基体 应 力 水 平 ,说 明 在 三 点 弯曲 载荷 
下 轴 向 纤维 棒 承 担 了 主要 载荷 ,基体 的 贡献 很 小 。 


图 3 三 维 编织 C/C 复合 材料 三 点 弯曲 染 宏 观 应 力 
Fig.3  Macro-stress of three-point bending beam of 


three-dimensional braided C/C composite 


+2.456e+07 


m ay 


(b) 下 表面 
图 4 三 维 编织 C/C 复合 材料 三 点 弯曲 梁 的 细 观 应 力 


Fig.4  Micro-stress of three-point bending beam of three-dimensional braided C/C composite 


势 基 本 吻合 。 在 加 载 初始 阶段 ,弯曲 载荷 - 挠 度 曲线 
基本 呈 线 性 。 当 载荷 接近 最 大 值 时 ,曲线 表现 出 轻 
微 的 非 线 性 ,这 是 由 基体 损伤 所 引起 的 。 随 后 , 载 
和 荷 - 挠 度 曲线 迅速 下 降 , 曲线 表现 出 弹 脆性 特点 ,出 
现 锯齿 状 波动 。 这 是 由 于 梁 从 上 下 表面 到 中 性 面 逐 


渐 破 坏 所 引起 的 载荷 重新 分 配 所 导致 。 选 取 弯 曲 载 
答 - 挠 度 曲线 上 载荷 最 大 值 作 为 梁 的 弯曲 失效 载荷 ， 
根据 式 (18 ) 可 以 计算 得 到 三 点 弯曲 梁 的 弯曲 破坏 
强度 o, 如 图 6 所 示 。3 个 试 样 的 测试 值 分 别 为 
129. 55 MPa 129.6 MPa 和 140. 9 MPa, ,测试 平均 值 
为 133. 35 MPa , 模型 预测 值 为 146. 85 MPa, 预测 值 
大 于 实验 值 。 这 可 能 是 由 于 数值 模型 中 没有 考虑 材 
料 制备 中 产生 的 空 际 、 微 裂纹 等 缺陷 ,以 及 纤维 束 / 
基体 界面 脱 粘 的 影响 。 
3pl 
2bh? 
式 中 :p HRAS Hl lUi LL 为 跨度 ;Lb WA AE; h 


(18) 


lh, gens 
pn CT 
> m lomi: 
qu b ; Tw 
©) oon 
M n 
co : im 
Ung. 
e 一 * 一 实验 -1 
le 一 实验 - 
" 一 0 一 实验 -3 

CD ENT 
eo 
co 0 0.5 1.0 1.5 2.0 2:5 
CN 挠 度 /mm 
Eds 三 维 纺织 C/C 复合 材料 当 载 和- 抄 度 曲线 
ae Fig. 5 Load-deflection curve of three-dimensional 


braided C/C composite beam 


140 90 146.85 


129.55 129.60 


chinaXiv 


实验 -1 实验 -2 实验 -3 预测 


q6 三 维 编织 C/C 复合 材料 弯曲 强度 实验 值 和 预测 值 


Fig.6 Experimental and predicted values of bending strength 


of three-dimensional braided C/C composites 
2.3 ”渐进 损伤 过 程 


图 7 为 利用 多 尺度 模型 模拟 获得 图 5 三 点 弯曲 
梁 载 荷 - 挠 度 曲线 上 工 、 工 、 亚 点 处 的 宏观 损伤 状态 。 
图 中 ，SDV1 ”表示 宏观 损伤 变量 ,是 含 损伤 的 单 胞 
均匀 化 刚度 矩阵 中 刚度 分 量 CS, Cd) 与 未 损伤 的 单 
胞 均匀 化 刚度 矩阵 中 刚度 分 量 C3 的 比值 , 含 损 伤 的 
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单 胞 均匀 化 刚度 C" (d) 是 根据 纤维 束 和 基体 的 损伤 
变量 d, ,利用 式 (9) 重 新 计算 后 获得 。“SDV1 值 在 
0 ~1 之 间 ,0 表示 没有 损伤 ,1 表示 完全 损伤 ,由 图 
可 见 , 在 三 点 弯曲 载荷 作用 下 , 梁 首 先 在 上 表面 中 间 
部 位 发 生 损 伤 , 随 着 载荷 的 增加 , 粱 上 下 表面 均 发 生 
损伤 ,损伤 的 区 域 逐 渐 扩 大 , 且 损伤 区 域 逐 渐 向 染 的 
中 性 面 扩展 ,上 下 表面 损伤 区 域 逐 步 贯穿 梁 的 厚度 
方向 截面 ,最 终 导致 梁 的 弯曲 失效 。 


SDV1 
(Avg: 100%) 


(b) IL zx 


+1.000e+00 
+9.000e-01 


+8.000e-01 
| +7.000e-01 
+6.000e-01 


+4.000e-01 
+3.000e-01 
+2.000e-01 
+1.000e-01 


+0.000e+00 

-7.660e-02 
7 
2J 

x? (c) TII 


图 7 三 维 编织 C/C 复合 材料 梁 三 点 弯曲 宏观 损伤 过 程 
Fig.7 Macroscopic damage process of three-point bending 
of 3D braided C/C composite beam 


为 了 验证 渐进 损伤 过 程 的 模拟 结果 ,进行 了 三 
维 编织 C/C 复合 材料 三 点 弯曲 实验 。 弯 曲 试 样 为 
60 mm x 14 mm x6.8 mm 的 长 方 体 , 跨 距 为 50 mm, 
三 点 弯曲 实验 采用 MTS 材料 实验 机 ,将 弯曲 试 样 放 
置 在 两 个 支撑 点 之 间 , 用 一 个 移动 压 尖 施加 准 静 态 
载荷 。 实 验 在 室温 条 件 下 进行 ,采用 位 移 加 载 控制 ， 
加 载 速度 1 mmymin, 记 录 试 样 加 载 点 处 的 位 移 和 载 
荷 。 图 8 为 三 维 编织 C/C 复合 材料 弯曲 试 样 和 弯 
曲 实验 。 


图 8 三 维 编织 C/C 复合 材料 弯曲 实验 
Fig.8  Three-dimensional braided C/C composites 


bending experiment 


048 


图 9 为 三 维 编织 C/C 复合 材料 梁 在 三 点 索 曲 
载荷 下 的 破坏 形 貌 。 由 图 9 可 见 ,三 维 编织 C/C 复 
合 材料 弯曲 破坏 的 部 位 和 区 域 大 小 与 数值 模拟 的 结 
果 符合 良好 。 


图 9 三 维 编织 C/C 复合 材料 梁 三 点 普 昌 破坏 形 肛 


Fig.9 Failure morphology of 3D braided C/C composite 


beam under three-point bending 


3 结 it 


CABE SERIE. 
.二 1) 基 于 多 尺度 渐进 展开 理论 ,结合 细 观 损伤 模 
型 建立 了 三 维 编织 C/C 复合 材料 结构 失效 的 宏 细 
尺度 分 析 方法 ,基于 商业 有 限 元 软件 ABAQUS 
JEBUT. UMAT 的 二 次 开发 ,实现 了 宏观 与 细 观 尺度 
之 闻 的 交互 信息 传递 以 及 多 尺度 损 损伤 模拟 。 
EASRA TRAR C/C 复合 材料 
失效 过 程 进行 了 模拟 ,预测 了 三 维 编织 C/C 复 
ms pu ul e E 
PORI, SUETA 13 0c EARE, WE 
ITERA Hou) Re UE AFB 
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